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We use the Kharzeev-Levin-Nardi (KLN) model of the low x gluon distributions to fit recent 
HERA data on Fl and F£ {F$). Having checked that this model gives a good description of the 
data, we use it to predict Fl and F% to be measured in a future electron-ion collider. The results are 
similar to those obtained with the de Florian-Sassot and Eskola-Paukkunen-Salgado nuclear gluon 
distributions. The conclusion of this exercise is that the KLN model, simple as it is, may still be 
used as an auxiliary tool to make estimates both for heavy ion and electron-ion collisions. 
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The small- x regime of QCD has been intensely in- 
vestigated in recent years (for recent reviews see, e.g. 

main prediction is a transition from the 
linear regime described by the DGLAP dynamics to a 
non-linear regime where parton recombination becomes 
important in the parton cascade and the evolution is gov- 
erned by a non-linear equation. At very small values of 
x we expect to observe the saturation of the growth of 
the gluon densities in hadrons and nuclei. One of the 
main topics of hadron physics to be explored in the new 
accelerators, such as the LHC and possibly the future 
electron-ion collider is the existence of this new compo- 
nent of the hadron wave function, denoted Color Glass 
Condensate (CGC) 0. 

The search for signatures of the CGC has been subject 
of an active research (for recent reviews see, e.g. (l|). Sat- 
uration models @, [H, 0, [E [E 0] can successfully describe 
HERA data in the small x and low Q 2 region. Moreover, 
some properties which appear naturally in the formalism 
of the color glass condensate have been observed experi- 
mentally. These include, for example, geometric scaling 
[1, [l(| and the supression of hig h pt h adron y ields at 
forward rapidities in dAu collisions [Til. 1T3, [l3l \l4 . Il5l. . 
However, it has been shown that both geometric scaling 
[ItJ and high px supression [l9[ can be understood 
with other explanations, not based on saturation physics. 

In view of these (and others) results we may conclude 
that there is some evidence for saturation at HERA and 
RHIC. However, more definite conclusions are not yet 
possible. In order to discriminate between these different 
models and test the CGC physics, it would be very impor- 
tant to consider an alternative search. To this purpose, 
the future electron-nucleus colliders offers a promising 
opportunity 0, EH, H [H, Hi] . 

The color glass condensate is important in itself as a 
new state of matter. However, apart from that, we need 



to know very well its properties since the CGC forms 
the initial state of the fluid created in nucleus- nucleus 
collisions. Any detailed simulation of a heavy ion col- 
lision needs a realistic Ansatz for the initial conditions. 
This would correspond to knowing accurately the unin- 
tegrated gluon distribution in the projectile and in the 
target. These distributions will presumably be known in 
the future, with the help of the results of deep inelas- 
tic scattering (DIS) off nuclei. Even before experiment, 
one can try to calculate the gluon density in the initial 
state of heavy ion collisions by numerically solving the 
classical Yang-Mills equations, as done in (25|. This is 
however very time consuming. Meanwhile, in practical 
applications we need to use use models for these distri- 
butions. Since these models are used as input for heavy 
numerical calculations, they have to be simple. 

One simple approach to saturation physics was devel- 
oped by Kharzeev, Levin and Nardi (KLN) in a series of 
papers [26| . where a simple model for the unintegrated 
gluon distribution was proposed. It was used in many 
phenomenological applications. In particular it was very 
successful when applied to hydrodynamical simulations 
[13, Hi and In [13 it was shown that hydrody- 

namical simulations with the KLN model initial condi- 
tions were able to the describe centrality, rapidity and 
energy dependence of charged hadron multiplicities very 
well. Moreover these simulations could reproduce the 
transverse momentum spectra of charged pions and also 
the centrality dependence of the nuclear modification fac- 
tors. 

In view of the success of the KLN model as an in- 
put for numerical simulations, we think that it would be 
interesting to confront it with recent DIS data and, if 
necessary, change it in order to get a better agreement 
with these data. Of course, the KLN will remain a phe- 
nomenological model to be replaced by something more 
fundamental in the future. However, in a refined and still 
simple version, it may be a very useful tool, capturing the 
essential physics of gluon saturation and parameterizing 



•Electronic address: babi@if.usp.br // duraes@mackenzie.br // navarrattfiiprxfeert^jfepi^ifitbiekcKpeiliEiental information. 



2 



In this paper we apply the KLN model to deep inelastic 
scattering, the domain where parton distributions have 
to be tested. If it fails badly in reproducing DIS data, it 
must be discarded. Since the KLN model gives only an 
Ansatz for the gluon distribution but says nothing about 
quarks, it is difficult to use it to make predictions for the 
most well known DIS observables, such as the structure 
functions F2 . We must then look for quantities which are 
dominated by the gluon content of the proton. These are 
the charm (bottom) and longitudinal structure functions, 
F 2 C (F2) and F L respectively. 

Let us first discuss charm production and its contribu- 
tion to the structure function. In the last years, both the 
HI and ZEUS collaborations have measured the charm 
component F$ °f the structure function at small x and 
have found it to be a large (approximately 25%) frac- 
tion of the total [3(|. This is in sharp contrast to what 
is found at large x, where typically F^jFi « 0(1O -2 ). 
This behavior is directly related to the growth of the 
gluon distribution at small- a;. 

In order to estimate the charm contribution to the 
structure function we consider the formalism developed 
in [3TI ] where the charm quark is treated as a heavy 
quark and its contribution is given by fixed-order per- 
turbation theory. This involves the computation of the 
boson-gluon fusion process j*g — * cc. A cc pair can be 
created by boson-gluon fusion when the squared invari- 
ant mass of the hadronic final state is W 2 > Am 2 . Since 

W 2 = Q ^~ x ^ +Mjy, where M N is the nucleon mass, the 
charm production can occur well below the Q 2 threshold, 
Q 2 w Am 2 , at small x. The charm contribution to the 
proton/nucleus structure function, in leading order (LO), 
is given by [13] 

^(x,Q 2 ,m 2 c ) = e 2 c ^fl f^CZtfrQgtof) , 
x 2vr J ax y y ' y 

(1) 

2 

where a = 1 + £ (£ = ^p-) and the renormalization scale 
ji is assumed to be either fi 2 — Am 2 or /1 2 = Am 2 + Q 2 . 
Cg 2 is the coefficient function given by 

C c g ,2 = { [2 2 + (1 - z? + 4 iz (1 - 3z) - 8 i 2 z 2 ] 
x lnH + 3[-l + 8z(l -z) -4£z(l - z)]} , (2) 

where 8=1- 4^ and H = ±±#. 

The dominant uncertainty in QCD calculations comes 
from the uncertainty in the charm quark mass. In this 
paper we assume m c — 1.2 GeV. In JT]) g(y,/i 2 ) is the 
glu on distribution, which is usually taken from the CTEQ 
j33j ] . MRST [IH or GRV ^ parameterizations. In what 
follows we shall use the KLN Ansatz: 

2 _/ - J w) s Q 2 ( l ~ X ) D « Q 2 <Q 2 , 

X9M) = \^SQl(l-x) D , Q*>Ql. (3) 

In the above expression S is the area of the target 
and a s is the running coupling a s = 12/[257r ln(Q 2 / A 2 )] 



(with A = 0.224 GeV). In © D = A and n is a con- 
stant parameter to be adjusted by requiring that the 
distribution xg{x,Q 2 ) satisfies the momentum sum rule 

Jq dx xg(x,Q 2 ) = p, where p is the value obtained with 
the GRV98 gluon density. Q s is the saturation scale given 
by: 

Q 2 s (x)=Ql(^f (4) 

where Q 2 Q = 0.34 GeV 2 , x = 3.0 x 10~ 3 and A = 0.25. 

In Fig. [T]we show F 2 as a function of x obtained with 
the above expression and compared to the ZEUS and 
HI data. Solid and dashed lines correspond to different 
choices for the renormalization scale. We can observe 
that there is a surprisingly good agreement between the 
KLN model and the data, especially considering that only 
minor changes in the parameters were made, with respect 
to those found previously in the analysis of RHIC data 
[2^] . We can also obtain a reasonable description of large 
Q 2 data, which is also surprising because the KLN model 
has no DGLAP evolution, being tuned to the low x and 
low Q 2 region of the phase space, where gluon satura- 
tion is expected to occur. In Fig. [2] we show the results 
for F$ and compare them with the HI data (36|. With 
the exception of the points with Q 2 — 200 GeV 2 , the 
agreement with data is similar to the one found for F£. 

New ex per imental HERA data on Fl have recently 
appeared [33] - The longitudinal structure function in 
deep inelastic scattering is one of the observables from 
which the gluon distribution can be unfolded. Longitu- 
dinal photons have zero helicity and can exist only vir- 
tually. In the quark model, helicity conservation of the 
electromagnetic vertex yields the Callan-Gross relation, 
Fl = 0, for scattering on quarks with spin 1/2 38]. This 
does not hold when the quarks acquire transverse mo- 
menta from QCD radiation. 

Instead, QCD yields the Altarelli-Martinelli equation 

M 

j? , ^ a s(Q 2 ) 2 f 1 dy r 8 2 
F L {x,Q ) = — £ x / —[-F 2 (y,Q ) 

+ 4£ e 2 (l--) 2 , 3 ( 2 ,,Q 2 )] , (5) 

q y 

cxpliciting the dependence of Fl on the strong coupling 
constant and the gluon density. At small x the second 
term with the gluon distribution is the dominant one. 
In Ref. [13] the authors have suggested that expres- 
sion §5§ can be reasonably approximated by Fl(x, Q 2 ) ~ 
0.3 ^-xg(2.5x, Q 2 ), which demonstrates the close rela- 
tion between the longitudinal structure function and the 
gluon distribution. 

In what follows we calculate Fl using the Altarelli- 
Martinelli equation, neglecting the F2 contribution and 
using ([3]) in ([5]). The results are presented in Fig. [3] 
where they are compared with the very recent HI data 
[37l ] and with the results obtained with the help of other 
standard gluon distribution functions. The agreement 
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between the KLN predictions and data is remarkably 
good. This is again surprising because the KLN dis- 
tribution is not expected to work so well at such large 
values of Q 2 . The good agreement indicates that the 
KLN distribution has a good asymptotic behavior and 
it is compatible both with the data and with the other 
standard gluon distributions. 

Having checked that the KLN distribution reproduces 
satisfactorily the existing DIS data on electron-proton 
collisions, we shall now use it to make predictions for 
electron-ion collisions. The expression for the nuclear 
charm structure function F£' is the same except for 
the change g(y, Q 2 ) -> g A {y, Q 2 ), where g A (y, Q 2 ) is ob- 
tained from Eq. with the replacements S — > S A = 
A 2 / 3 S and Q 2 -> Q 2 S A , where Q 2A = A^ 3 Q 2 S . In order 
to estimate the strength of non-linear effects in eA pro- 



cesses, we can compute the linear contribution to F^ us- 
ing only the second line of Eq. . We call this p^" mear ) 

and compare it with F^ \ where the latter is calcu- 
lated with both lines of Eq. © and thus includes non- 
linear effects. The ratio pc(f utl ) / pc( lmear ) as a function 
of x is shown in Fig. 0] for A — 208 and for several values 
of Q 2 . 

The deviation of this ratio from unity shows the im- 
portance of non-linear effects. As expected, for large x 
and for large Q 2 there are no saturation effects. In fact, 
saturation effects are only noticeable at x < 10 ~ 6 and for 
small values of Q 2 . This result confirms the findings of 
[4l| . where an estimate of saturation effects in eA colli- 
sions performed with the color dipole approach also led 
to the conclusion that they are only marginally visible in 
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FIG. 3: Longitudinal structure function Fl computed with 
the KLN model and with the gluon distributions taken from 
CTEQ6L (dash dotted line), MRST98 (dotted line) and 
GRV98 (dashed line). Data are from Ref. [13 



i 7 ^'" 4 . Having an idea of where saturation effects could be 



relevant, we can compute an observable quantity, which 
is R c (x, Q 2 ) = F^' A (x, Q 2 )/AF^' p (x, Q 2 ). The deviation 
from unity in this ratio is an indication of saturation 
physics. A depletion in this ratio is called "shadowing" , 
whereas an enhancement is called "anti-shadowing" . 
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In Fig. [5] on the top panel we calculate R c and esti- 
mate the magnitude of shadowing, which can be of 50 % 
at very low (but still reachable) values of x and Q 2 . The 
equivalent ratio for the longitudinal structure functions 
R L (x, Q 2 ) = Ff (x, Q 2 )/AFl(x, Q 2 ) is shown in the bot- 
tom panel of Fig. [5] for the same choices of Q 2 . We 
observe that significant non-linear effects start to appear 
at larger values of a;. Rl seems thus more promising than 
R c . Following (4l| we compare the normalized ratios R c 
and Rl obtained with the KNL model with the same 
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FIG. 6: Ratio R c calculated with the KLN, DS and EPS 
nuclear gluon distribution functions for A — 208 and Q 2 = 2 
(top) and Q 2 = 11 (bottom) GeV 2 . 
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ratios computed with the standard collinear factoriza- 
tion approach with nuclear parton distribution functions 
(nPDF's). We take two extreme cases, one with almost 
no shadowing at all, based on the nPDF's of de Florian 
and Sassot (called here DS) [13] and one with maximum 
shadowing based on the Eskola, Paukkunen and Salgado 
nPDF's (called here EPS) [H[. 

The ratios R c and are shown in Figs. [S] and [7] re- 
spectively. From the figures we see that the KLN model 
interpolates between the two extreme parameterizations, 
DS and EPS, being closer to the latter. This is expected, 
since the EPS gluon distribution comes from a fit of world 
data where BRAHMS data on forward particle produc- 
tion were included. Both KLN and EPS take RHIC data 
into account. 

In summary, we have used the KLN model for the low 
x gluon distributions, slightly changing the parameters 
fixed from previous analysis, to fit HERA data on Fl 
and i<2- Having checked that this model gives a good 
description of the data, we have used it to predict Fl and 
F% to be measured in electron-ion collisions. The results 
are close to those obtained with the DS and EPS nuclear 
gluon distributions. The conclusion of this exercise is 
that the KLN model, simple as it is, may still be used as 
an auxiliary tool to make estimates both for heavy ion 
and electron-ion collisions. 
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